Linear local modes induced by intrinsic localized modes in a monatomic chain 
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A theory is developed to describe the effect of an intrinsic localized mode (ILM) on small vibrations 
in a monatomic chain with hard quartic anharmonicity. One prediction is the appearance in the 
chain of linear local modes nearby the ILM. To check this result, MD calculations of vibrations 
under strong local excitation are carried through with high precision. The results fully confirm the 
prediction. 
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I. INTRODUCTION 

The realization [iHl] that there can exist stable 
strongly localized excitations in perfect classical anhar- 
monic lattices has led to theoretical studies exploring a 
variety of possibilities [341^ . The frequency of such a lo- 
calized vibration depends on its amplitude and lies out- 
side the phonon spectrum. Such excitations are called 
intrinsic localized modes (ILMs) 0, to emphasize their 
similarity in appearance to defect impurity modes, or 
called discrete breathers or discrete solitons [H, [13] 
to make a connection to solitons in continuous systems. 
Some experimental evidence for the existence of intrinsic 
localized modes in 1-D lattices in microsco pic and macro- 
scopic lattices has been demonstrated |15l - [23 |. 

It is to be expected that the appearance of an ILM 
would change the local properties of the lattice includ- 
ing the local phonon dynamics. This back reaction on 
the phonon spectrum should have physical consequences 
since the ILM could induce local modes outside the plane 
wave spectrum. The presence of these additional res- 
onances and their dependence on the amplitude of the 
ILM should add some complexity to the energy relax- 
ation rate of the ILM to the phonon bath. In a different 
direction, since an ILM can move through the lattice it is 
expected that such trapped local modes would also move 
or, more likely, tend to inhibit the translational motion 
of such localized energy. 

In this communication we examine analytically the 
small amplitude vibrations of a 1-D nonlinear chain with 
intersite coupling in the presence of an ILM. The results 
show that an ILM with sufficient nonlinear amplitude 
stabilizes the appearance of linear local modes (LLMs) 
above the top of the plane wave spectrum. Next molec- 
ular dynamics (MD) simulations are used to verify these 



findings. First an ILM is formed and the power spectrum 
calculated, next the initial conditions are changed and 
very small amplitude shifts associated with a LLM are 
added to the neighboring atoms and the power spectrum 
calculated again. The cases when the new spectra contain 
one additional weak peak above the phonon spectrum are 
studied in detail. The frequencies of the additional peaks 
do not depend on their amplitudes but they do depend 
on the ILM amplitude consistent with the linear require- 
ment. Good agreement is obtained between the positions 
of the additional peaks and the theoretically calculated 
frequencies of the LLMs. 



II. THEORY 

Let us consider an anharmonic monatomic chain with 
interactions of nearest-neighboring atoms. Taking into 
account only quartic anharmonicity, the equation of mo- 
tion of atoms in this model reads 



E 

/=n±l 



[fc2 (C/n' - C/n) + fc4 (C/n' 



(1) 
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where C/„ are the reduced displacements of atoms located 
at the site n of the chain, the subscripts n indicate the 
number of the site, /c2 and are the parameters of har- 
monic and anharmonic springs, whereat k2 = de- 
termines the top phonon frequency LOm- If A:4 > then 
ILMs may be excited in the chain with the frequencies 
above the allowed phonon spectrum 0, 0] . 

To describe small vibrations of the chain in presence 
of an ILM, we add to [/„ an infinitesimal displacement 
g„: U'^{t) — Un{t) + qn{t)- For the displacements of 
the ILM we take C/„(t) — An cos ojLt, where An is the 
amplitude parameter of the ILM which may be both, 
positive or negative (the small contribution of higher- 
order harmonics is omitted). The shift C/^(t) also satisfies 
Eq. ID) . Subtracting from this equation the equation for 
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Un{t) we get 

9n = ^ [k2+3k4{An-An'fcOS^UJLt]{qn'-qn)- (2) 
n'— nil 

One can divide Qn into two parts: 1) the shifts Qo^n 
describing small variations of the ILM (they have been 
considered in j23l, [24] ) 2) all other shifts qi.n- Here we 
are interested in the stable solutions of equation ([2]) which 
correspond to the latter modes. These modes are orthog- 
onal to (i.e. independent from) the ILM. We can take the 
orthogonality condition into account if we add in Eq. ^ 
the term A An An' Qn' , where the Lagrange multiplier 
A should be found from the orthogonality condition of 
the LLM under consideration. However, if one consid- 
ers the linear modes of symmetry different from the ILM 
then there is no need to add the term oc A, since the 
orthogonality conditions are already fulfilled. 

To identify LLMs, we first apply the rotating wave 
approximation which corresponds to the replacing in 
equation ([2]) the factor cos^ uj^t = (1 -|- cos2wit)/2 
by 1/2 (i.e. we neglect the oscillating in time term 
(1/2) cos2a;i,t). In this approximation the equation of 
motion of the linear modes takes the form: 



^wn' jQn' 1 



(3) 



where Dnn' — k2{2Snn' —^n±i,n') IS the dynamical matrix 
of the perfect monatomic chain. 



[Snn' {{An ~ Ai+l) 



-^n±l,n 



{An~An'f\ + 



\An An' 



0^) 



(4) 



is the perturbation of the dynamical matrix. 

The effect of the perturbation can be found by the Lif- 
shitz method [l^ . Thus a linear local mode / exists if the 
imaginary part of the Green's function of the perturbed 
lattice has a pole at cj; outside the allowed phonon spec- 
trum. The latter function can be found from the equation 

G(cj) ^ (/-G(")(w)w)"'G(°n^), (5) 
where G, f , and are matrices; v is given by Eq. 

Gi°i-(^) = hp(a.)]l"-"'l/c^v/^^, (6) 

is the (n, n')-component of the Green's function matrix 
of the perfect chain [2^, p(bj) = (w— — 1)^ < 1 (the 
units uJm = 1 are used). The amplitude parameters of a 
LLM are given by the relation 



an, I = G„„i(aj/)/G„i„i(aj;). 



(7) 



The time oscillatory terms in Eq. ([2]), neglected in 
the rotating wave approximation, can lead to new effects 
not found in the vibrations of a lattice with a static de- 
fect. This difference follows from the Floquet theorem 



according to which in a periodically time-dependent sys- 
tem with the period T, besides the excitations with the 
frequency ioi, there exist also excitations with the fre- 
quencies I ± w; -|- 2TrN/T\, where N is an integer. In our 
case T — tt/lul] the repetitions of the frequency result 
from the time-dependent part of the anharmonic inter- 
action of the LLMs with ILMs. 

The time oscillatory terms in Eq. ^ can also cause 
a renormalization of the frequencies of linear modes. 
This follows from the fact that these terms oscillate 
in time with the frequencies 2ujl—u! and 2ujl+^^ (here 
uj is the frequency of a small vibration), one of which, 
2lul — i-^i, may be comparable with the perturbation v. 
As a result, the frequencies of some linear modes may 
acquire complex values, which will cause these modes to 
become unstable (27n3Q] . 

To find whether a LLM is stable or not we apply in 
Eq. the expansion g„ = J2j ^njXj, where e„j is the 
normalized contribution of the atom n to the normal co- 
ordinate Xj of the perturbed lattice. We get 



Xi = Ujfxj + C0s2uJLt ^ ^ e„; Wnn'en'j Xj. 



(8) 



The effect of the time oscillatory term in Eq. is most 
important for the modes with frequencies close to w^. 
Supposing that only the LLM under consideration has 
such a frequency, we can neglect in the right-hand-side 
of this equation the terms j ^ I. In this approximation 



xi = t^f (1 + hcos2ujLt) xi. 



(9) 



where h — 



I Wn 



Cn'i- (Here the inhomoge- 



neous term oc XL{t) cos 2ujLt also has been neglected; this 
term describes the infinitesimal forced oscillations with 
the frequencies lol and SojLresulting in the infinitesimal 
shift of the ILM phase.) Equation © is the Mathieu 
equation describing the parametric resonance [sij . The 
term cx h may result in a) renormalization of the fre- 
quency of the LLM, or b) in the modes instability; in the 
latter case the renormalized frequency is complex. Be- 
low we consider the renormalization and the instability 
conditions numerically. 



A. LLMs induced by even ILM 

Both even and odd ILMs can exist in the nonlinear 
monatomic chain [3]. Here we consider an even ILM to 
be the effective defect in this lattice calculation, since this 
mode is stable with respect to small translational fluctu- 
This choice will permit calculations of ILMs 



ations 

+ LLMs with high precision. The amplitude patterns 
for such an ILM for different values of the dimensionless 
nonlinear parameter, ^4^0/^2 , are given in rows 3 and 
4 of Table 1. 

An even ILM can induce the appearance of both, odd 
and even LLMs. In the case of odd LLMs the A- mul- 
tiplier in Eq. ^ equals zero. Then, if one considers 
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an ILM localized of six neighboring central atoms, the 
perturbation matrix v (see Eq. (|3])) is the six-range sym- 
metric matrix with the following nonzero elements: 

voo = vii = 70 + 7i, Vol = -70, V12 = -71; 

V22 = 71 +72, «23 = -72, V33 = 73, (10) 

and with few matrix elements with negative index sat- 
isfying the relations v^n-n' = Vn+in'+i {n > 0) . Here 
7„ « 3fc4(A„ — A„+i)^/2fc2 are the changes of the springs 
(in dimensionless units) . The effect of an even ILM upon 
even phonons is described by the perturbation matrix 
given by Eq. ([4]), which is equal to the previous per- 
turbation matrix v plus the term XAnA^i ; the value of 
A must be found self-consistently from the orthogonality 

condition Yl,n (in,evnAn = 0. 

Inserting the matrix v into Eqs. ([5]) and (O gives the 
frequencies and relative amplitudes of the LLMs. The 
parametric resonance parameter h, which determines the 
frequency correction equals 



^(|e„| + |e„+i|)^(|A„| + K+i|)'. (11) 



,3X4^1 

* n>l 

Table 1. Analytical calculations of odd and even 
LLMs for a lattice containing an even ILM with a range 
of nonlinear parameter ^4^0/^2 values. Given are the 
relative ILM frequency oJi/wm, the ILM amplitudes of 
the four next-to central particles with = — y4„+i, 
the relative frequency of the odd and even LLM ujodd/'-^m 
and ojevn/^m, the amplitudes of these modes a^n.odd = 

On+l.odd and a^n.evn — ^0'n+l,evn {n > 0); Ai — Aq — 1; 

01, odd — ai^evn — 1. (The frequencies of the odd LLMs 
are corrected for the time oscillatory perturbation; h and 
her are the corresponding parameters of this perturba- 
tion.) 



k4Al/k2 


0.50 


0.75 


1.00 


1.25 


1.50 


1.75 


2.00 




1.36 


1.51 


1.64 


1.77 


1.88 


2.00 


2.10 


-A2 


0.199 


0.170 


0.153 


0.141 


0.133 


0.127 


0.121 


A3 


0.042 


0.027 


0.019 


0.015 


0.012 


0.010 


0.008 


^oddl^m 


1.288 


1.418 


1.486 


1.58 


1.677 


1.78 


1.86 


Cil,odd 


0.729 


0.784 


0.819 


0.844 


0.862 


0.876 


0.887 


0-S,odd 


0.347 


0.262 


0.212 


0.179 


0.156 


0.138 


0.124 


hodd 


0.456 


0.560 


0.629 


0.679 


0.711 


0.746 


0.770 


hodd,cr 


0.52 


0.63 


0.77 


0.84 


0.89 


0.94 


0.98 


^evn 1 


1.009 


1.031 


1.063 


1.102 


1.142 


1.179 


1.219 




0.383 


0.308 


0.273 


0.252 


0.238 


0.228 


0.220 


0-3, evn 


1.143 


0.822 


0.639 


0.522 


0.442 


0.384 


0.340 


A -102 


0.70 


0.83 


0.87 


0.92 


0.96 


0.98 


1.01 



In Table 1 the calculated values for odd and even LLMs 
are presented for different values of the nonlinear parame- 
ter of the ILM, A:4Ag/fc2. Comparing the LLM amplitude 
patterns with the amplitudes of the ILM shows that they 
are mutually orthogonal and the ILM and the LLM be- 
long to the different degrees of freedom. Note that the 
frequency of the odd LLM is rather close to ojl. For 
this mode consideration of the time oscillatory term in 
Eq. ([9]) is important. The numerical calculations of the 
Mathieu equation ([9|) shows that in all cases the paramet- 
ric resonance parameter h is less than its critical value 
her for the instability. The frequency of the even LLM 
is significantly different from wz, so that the correction 
caused by the time oscillatory term is now small. 



III. MD SIMULATIONS 

To perform MD simulations of ILMs and LLMs, we 
integrate numerically the equation of motion given by 
Eq. We calculate the vibrational amplitudes and 

velocities of atoms of the finite chain with 10^ atoms in 
10^ time points. The time interval includes ~ 10'^ periods 
of the ILMs. We also calculate the vibrational power 
spectrum. To decrease the background in the spectrum, 
Fourier transformations are carried out between the time 
points corresponding to the same phase of vibrations. 
For the initial condition, we used the ILM displacements 
of the six central atoms from their equilibrium positions. 
The amplitude difference of the two central atoms is fixed 
at unity with k2 = 100 {uj,n = 20). In this case the actual 
values of k4 for ILMs with frequency < 2L0m are of the 
order of k2 ■ The initial shifts of other four central atoms 
are selected so that only the ILM is excited. The resulting 
amplitude parameters An of six central atoms are those 
shown in Table 1. 

Examples of the calculated power spectrum of an even 
ILM are given in the bottom panels of Fig. 1. Two 
different values of the nonlinear parameter are shown. 
Each spectrum has a single peak with a frequency 
i^L > Not shown are the much weaker peaks asso- 
ciated with the higher odd harmonics at Suj^, 5wi, etc. 
Hence the bottom panels demonstrate that our initial 
conditions correspond solely to the excitation of an ILM 
(at least with the accuracy 10"^''). 

When small additional amplitude shifts appropriate to 
the odd or even symmetry LLM presented in Table 1 
are added to the ILM amplitude pattern, then additional 
small peaks appear in the power spectrum, one 
at w < Wi, and an even smaller symmetrically situated 
peak at a; > lol. The middle panels in Fig. 1 show the 
case for an odd LLM, with its amplitude pattern, and 
the top panel for an even LLM, again with its amplitude 
pattern. One finds analogous sideband spectra for other 
A:4v4g/A:2 values. These data show that the frequencies 
of the new spectral peaks depend on the amplitude and 
hence the frequency of the ILM. At the same time the 
solid and dotted traces presented in Fig. 1, where the 
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FIG. 1: Power spectra of ILMs for two different nonlinear pa- 
rameter values with and without LLMs. Left column of fig- 
ures: the main spectral feature is an even ILM with k^^A^jk^ 
= 0.75, right column of figures: same even ILM but with 
fc4^o/fe = 1-5. Bottom panels: the unperturbed even ILM; 
middle panels: an ILM perturbed by an odd LLM; top panels: 
an ILM perturbed by an even LLM. Although the dotted and 
solid lines differ by lO'^ in the LLM amplitude its frequency 
remains unchanged. The amplitude patterns of the ILMs and 
LLMs are shown by arrows. 



LLM amplitude is varied by 10^, demonstrate that these 
LLM frequencies do not depend on their own amplitudes 
as expected for a linear response (the amplitudes of the 
LLMs are ~ 10^ and ~ 10^ times less than the amplitude 
of the ILM). The = 2J7i — w; of the symmetrically 
situated weak peak indicates that this spectral feature is 
the four wave mixing response of the LLM with the ILM. 



IV. DISCUSSION AND CONCLUSIONS 

There is value in a quantitative comparison of 
the analytical results and MD simulations for the 
LLMs. Figure 2 presents the frequencies of the LLMs ob- 
tained by both these methods. Inspection of these results 
shows that there is excellent agreement between the two 
methods. Note that the time oscillatory terms provide a 
significant correction for odd LLM case but are insignif- 
icant for the even LLM. The reason is in the relatively 
small difference of wl — Wodd as compared to — (jJevn ■ 

We have presented a theory which allows one to de- 
scribe the effect of an ILM on phonons in a nonlin- 
ear monatomic chain. The prediction of linear localized 
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FIG. 2: The dependence of the frequencies of the odd and even 
LLMs on the dimensionless nonlinear parameter k/^A^jk^ of 
the ILM. The frequencies of the odd LLM are given with (dot- 
ted line) and without (dash-dotted line) taking into account 
the time-oscillatory term hcos2ijJit in Eq. ((9]). The ILM fre- 
quency dependence is also shown. 



modes above the top of the plane wave spectrum is in 
good agreement with MD simulations. Basically the ap- 
pearance of the ILM changes the nearby nonlinear spring 
constants sufficiently so that linear local modes can also 
appear. The resulting lattice perturbation produced by 
the ILM is analogous to that associated with a force con- 
stant defect in a linear lattice in that a small number 
of degrees of freedom can be strongly perturbed but all 
phonons are perturbed to some extent. Looking ahead 
to more realistic diatomic lattices involving two body po- 
tentials |9|] , we expect that a variety of LLM possibilities 
may appear. These could include modes of the local, gap, 
and resonant types as well as tunneling states, all made 
possible by the fundamental combination of nonlinear- 
ity and lattice discretness. One may anticipate that, in 
analogy with the previously studied force constant defect 
cases, ILM- induced IR and Raman-activity may occur 
throughout the entire phonon spectrum [Sa fssj . 

There are also some very interesting differences be- 
tween the properties of extrinsic and intrinsic local- 
ized excitations evident from our study of this simple 
model system. One is that an ILM introduces a periodi- 
cally time-dependent perturbation that supports nonlin- 
ear mixing processes between the ILM and the LLMs and 
a second is that for a moving ILM the "effective defect 
space" with its associated LLMs will travel together with 
the nonlinear excitation. 
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